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' Abstract 

We consider the motion by mean curvature of an n-dimensional 

I grapli over a time-dependent domain in R", intersecting R" at a con- 

^! ■ stant angle. In the general case, we prove local existence for the cor- 

I responding quasilinear parabolic equation with a free boundary, and 

■ derive a continuation criterion based on the second fundamental form. 
I If the initial graph is concave, we show this is preserved, and that the 

I— !■ solution exists only for finite time. This corresponds to a symmetric 

Oh I version of mean curvature motion of a network of hypersurfaces with 

■ triple junctions, with constant contact angle at the junctions. 

j3 ' !• Time-dependent graphs with a contact angle condition. 

S '. 

I— li We consider a moving hypersurface T,t in M"^-'^, with normal velocity 

^ I equal to its mean curvature, assumed to be a graph over a time-dependent 

^ ■ open set D{t) C M" (not necessarily bounded, or connected.) The (properly 

CN , embedded) {n — l)-submanifold of intersection: 

ON . 

T{t) = ^tnW = dD{t) 

ly-^ . is a 'moving boundary'. Along T{t) we impose a constant-angle condition: 

{N,en+i)\r{t) = /3, 

o ■ 

•• . where < /3 < 1 is a constant and N is the upward unit normal of S^. 

. ^ I 'Mean curvature motion' is defined by the law: 

^ ■ Vn = H, 

where V/v = {V, N), with V = dfF the velocity vector in a given parametriza- 
tion F{t) of Tit {V depends on the parametrization, while V/v does not). A 
particular parametrization yields 'mean curvature flow': 

dtF = HN. 

For graphs, it is natural to consider 'graph mean curvature motion': if = 
graph w{t) for a function w{t) : D{t) M, imposing {dtF,N) = H with 
F{y,t) = [y,w{y,t)] for y G D{t), we find: 

w^ = y^l + \Dw\'^H 
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(and the velocity is vertical, dtF = wten+i)- With the contact angle condi- 
tion, we obtain a free boundary problem for a quasilinear PDE: 



2 



f wt = g'^ {Dw )wij in D {t), 

\ w = Q, fj^/l + \Dw\^ = 1 on dD{t), 

where g^^{Dw) = 5^^ — WiWj/(l + |Z)ti;p) is the inverse metric matrix. 

Remark 1.1. It is easy to see that the const ant- angle boundary condition 
is incompatible with mean curvature flow parametrized over a fixed domain 
Dq: on dDfy we would have (F, en+i) = 0, leading to {dtF,en+i) = 0, 
incompatible with dtF = HN and (AT, e„+i) = p. If we parametrize over 
time-dependent domains, mean curvature flow and graph m.c.m. lead to 
identical normal velocities for the moving boundary (sec section 2.) 

To establish short-time existence (in parabolic Holder spaces) we will 
work with a third realization of the motion, defined over a fixed domain: 

F{t) : Do M"+^ F{x, t) = [ip{x, t), u{x, t)] e x M, 

where (p{t) : Dq D{t) is a diffeomorphism and F is a solution of the 
parabolic system: 

Ft = g'={DF)F,j, 

where gij = 6ij + {Fi,Fj) is the induced metric on and g^^ is the inverse 
metric matrix. 

In the first part of the paper (sections 3 to 8) we prove the following 
short-time existence theorem (on Q := Dq x [0, T]: 

Theorem 1.1. Let Sq C M"+^ be a (7^+° graph over Dq C W satisfying 
the contact and angle conditions at dDo. There exists a parametrization 

-^0 = S C'^^°'{Do) of Eq, T > depending only on Fq and a unique 

solution F G C2+"'i+"/2(Q^; M"+i) of the system: 

r dtF-g'^{DF)didjF = 0, 
\u\aD,=0, iV"+i {D^, Du) laoo = /3, 

with initial data Fq and satisfying, in addition, the 'orthogonality conditions' 
at ODq (described in section 3.) 

The system and boundary conditions are discussed in more detail in 
section 3. Sections 4, 5, and 6 deal with compatibility at t = 0, linearization 
and the verification that the boundary conditions satisfy 'complementarity'. 
In particular, adjusting the initial diffeomorphism (^0 to ensure compatibility 
(section 4) leads to the 'loss of differentiability' seen in theorem 1.1. The 
required estimates in Holder spaces for the linearized system are described 
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in section 7, and the proof concluded (by a fixed-point argument) in section 
8. While the general scheme is standard, due to the particular boundary 
conditions adopted many details had to be worked out from first principles. 
Free boundary-type problems for mean curvature motion of graphs have 
apparently not previously been considered. 

We describe the evolution equations in the rotationally symmetric case 
in section 9 (including a stationary example for the exterior problem) and 
the extension to the case of a graph motion intersecting fixed support 
hypersurfaces orthogonally in section 10. 

The original motivation for this work was to establish (by classical parabolic 
PDE methods) existence-uniqueness for mean curvature motion of networks 
of surfaces meeting along triple junctions with constant-angle conditions. 
One can use a motion of graphs with constant contact angle to produce 
examples of 'triple junction motion': three hypersurfaces moving by mean 
curvature meeting along an (n — l)-dimensional submanifold $](t) so that 
the three normals make constant angles (say, 120 degrees) along T(t). The 
simplest way to do this is by reflection on M", so the hypersurfaces are T,t, 
and R" — D{t). If = graph w^t) with w > 0, the system is embedded in 
j^n+i^ This is mean curvature motion of a 'symmetric triple junction of 
graphs'. 

Short-time existence holds for general triple junctions of graphs moving 
by mean curvature with constant 120 degree angles at the junction, provided 
a compatibility condition holds along the junction (see section 15). The idea 
of proof is similar to the one given here; since the details are easier to un- 
derstand in the symmetric case, we decided to do this first. In addition, in 
the present case it is possible to go a lot further towards a geometric global 
existence result. Motivated by recent work on 'lens-type' curve networks [5], 
in the second part of the paper (sections 11-14) we consider continuation cri- 
teria and preservation of concavity. Since we chose to develop these results 
for graph motion with a free boundary, although the general lines of proof 
(via maximum principles) have precedents, many details had to be devel- 
oped anew. For example, section 13 contains an extension of the maximum 
principle for symmetric tensors with Neumann-type boundary conditions 
given in [9], which in our setting allows one to show preservation of weak 
concavity in general. The results obtained in sections 11-14 are summarized 
in the following theorem. 
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Theorem 1.2. Let Tmax be the maximal existence time for the evolu- 
tion. Assuming Tmax < oo, the second fundamental form h is unbounded at 
the junction Tt, as t ^ Tmax- 

limsup(sup = oo. 

If the mean curvature of the initial hypersurface is strictly negative (sup^^^ H = 
Hq < 0), then Tmax is finite. If Sq is weakly concave (/i < at t = 0), this 
is preserved by the evolution. 

The expected global existence result is that, assuming weak concavity, 
diam{T,t) as t ^ Tmax- 

Acknowledgments. It is a pleasure to thank Nicholas Alikakos for origi- 
nally proposing to consider the problem of mean curvature flow for networks 
of surfaces meeting at constant angles, and for his interest in this work. 
Most of the work on short-time existence was undertaken during a stay at 
the Max-Planck Institute for Gravitational Physics in Golm (January-June, 
2007); I am grateful to the Max-Planck Society for supporting the visit, 
and to Gerhard Huisken, director of the Geometric Analysis group, for the 
invitation. Finally, thanks to Mariel Saez for communicating the results of 
the Lens Seminar ([5]) and of her recent work on mean curvature flow of 
networks (partly in collaboration with Rafe Mazzeo, f8]). 

2. Normal velocity of the moving boundary. The evolution is nat- 
urally supplied with initial data Sq, a graph meeting M"'"'"^ at the prescribed 
angle. Since we are interested in classical solutions in the parabolic Holder 
space we expect an additional compatibility condition at t = 0. 

We discuss this first for graph m.c.m. w{y,t). 

Denote by T(t) a global parametrization of dD{t) (with domain in a 
fixed manifold, and 'space variables' left implicit). Differentiating in t the 
'contact condition' 'w(T{t),t) = 0, we find: 

wt + {Dw,t{t)) =0. 

Denote by nt the unit normal vector field to T{t), chosen so that {ut, Dw) > 
0. The contact condition also implies the gradient of w is purely normal: 

Dw\dD{t) = {Dntw)nt. 

Combining this with the angle condition, and bearing in mind that -Dnt 'U^|r(t) > 
0, we find: 

Dn,w = ^ on dD{t), po := ^1-/32. 
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Thus, on dD{t): 

= Vl + {D^,wYH = wt = -{t{t),nt)Dn,w = -i'n{t)j, 

and we find the normal velocity of the moving boundary (independent of 
the parametrization of Fj): 

which in particular must hold at i = 0. Note that we don't get a 'com- 
patibility condition' in the usual sense (of a constraint on the 2-jet of the 
initial data), but instead an equation of motion for the moving boundary. 
(Later, in the fixed-domain formulation, we will have to deal with a real 
compatibility condition). 

Now consider mean curvature flow parametrized over a time-dependent 
domain ^{t), with the boundary conditions: 

{F{x, t), en+i) = 0, {N, en+i) = P, xe &D{t). 

Let u = ut he the inner unit normal to P(t). Suppose S{9,t), 6 G S*"^^, 
parametrizes dD{t)] thus the 'junction' 9Sf is parametrized by V{9,t) = 
F{t,S{9,t)), and (denoting partial t derivatives with a dot): 

t{e, t) = dtF{t, S{e, t)) + dF[S{9, t)] = HN + (<S • z/)5^F 

(where we used the fact that Ot-F = for any r G T&D{t).) Thus, using 
{N,n) = -po: 

tn := t{9,t)-n = -poH+iS ■ v){d,F,n). 

On the other hand, from {F{t, S{9, i)), e„+i) = 0, we find by differentiation: 

Hf3+{S-v){d,F,en+i) = ^, 

OT S ■ V = —Hj3/{dj^F,en+i)- Letting T := -j^^ (tangent to S( at the 
interface), we have: 

f„ = -g(/3o + /3 ,^^'"\ ). 
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Denoting by A^' = — Pen+i the component of A^, we clearly have 
n = -{l/Po)N', so (T,n) = -(l//3o)(T, iV') = (/3//3o)(r, e„+i), and we 
conclude: 

/?2 1 

Po Po 

as before. 

Remark 2.1. This is not unexpected, if we accept there is a reparametriza- 
tion connecting the two motions, respecting the boundary conditions. That 
is, the ODE argument in [Ij should also work in the presence of boundary 
conditions and moving boundaries. 

Remark 2.2. We remark that for more general (non-symmetric, non-flat) 
triple junctions with 120 degree angles, the condition: 

+ H'^ = on T{t) 

must hold at the junction (for graphs, oriented by the upward normal), which 
in particular gives a geometric constraint on the initial data, for classical 
evolution in C^+^'i+^Z^. This is automatic in the symmetric case {w = 
—w^), since //^ = and = trgicPw^ for 7=1,2. 

3. Choice of 'gauge'. It is traditional in moving boundary problems 
to parametrize the time-dependent domain D(t) of the unknown w{y,t) by 
a time-dependent diffeomorphism: 

y = ip{x,t), ip{t):Do^D{t), 

and then derive the equation satisfied by the coordinate-changed function 
from the equation for w (see e.g. [6] or |10j). Motivated by the work on 
curve networks ([7]) we will, instead, consider a general parametrization: 

F -.DqX [0, T] M"+\ F{x, t) = [ip{x, t),u{x, t)] G M" x M 

and derive an equation for F directly from the definition of mean curvature 
motion: 

{dtF,N)=H. 

(We'll still assume ip{t) : Dq — > D{t) is a diffeomorphism.) The first and 
second fundamental forms are given by: 

g,j = {Fi ,Fj), A{Fi ,Fj) = {Fij , N) . 
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(Notation: DF = FiSi, D'^F{ei, ej) = Fij, (cj) is the standard basis of M"+^.) 
The mean curvature is the trace of A in the induced metric: 

H = {g'^{DF)Fi,,N). 

The equation for F is: 

{dtF-g'^{DF)Fij,N) =0. 

There is a natural 'gauge choice' yielding a quasilinear parabolic system: 

dtF-g'^{DF)Fi,=0. 

We will sometimes refer to this as the 'split gauge', since in terms of the 
components F = u] we have the essentially decoupled system: 

f dtu- g'^{D(f,Du)uij = 0, 
\ dt(p - g'^ (Dip, Du)ipij = 0. 

The splitting is useful to state the boundary conditions: 

u^dDo — ('contact'), 
iV"+i(I)^,Z)u)|aBo =/3 ('angle'). 

We immediately see there is a problem, since we have 2 scalar boundary 
conditions for n + 1 unknowns (and no moving boundary to help!) Our 
solution to this is to introduce n — 1 additional 'orthogonality conditions' at 
the boundary for the parametrization (p{t). We impose: 

{Dr<f,Dn(p)\dDo = 0, 

for any r E TODq, where n is the inward unit normal to Dq. 

Geometrically, the 'orthogonality' boundary condition has precedent in 
a method often adopted when dealing with the evolution of hypersurfaces 
in R"+^ intersecting a fixed n-dimensional 'support surface' orthogonally 
(see e.g. [H]): one replaces vanishing inner product of the unit normals 
(a single scalar condition) by a stronger Neumann-type condition for the 
parametrization, corresponding to n — 1 scalar conditions. (More details are 
given in Section 10.) 

The system must also be supplied with initial data. We assume given an 
initial hypersurface Sq, the graph of a C^"^" function uo{x) defined in the 
(j3+a domain Dq C M". (The reason for this choice of differentiability class 
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will be seen later.) It would seem natural to set cpo = Idoo, but this causes 
problems (related to compatibility; sec Section 4 below). We do require the 
1-jet of ipo at the boundary to be that of the identity: 

'PO\dDo = M Dipo^Qjj = I. 

(In particular, the orthogonality condition holds at i = 0.) 

We need a more explicit expression for the unit normal, and for that we 
use the 'vector product': 



N{D<f,Du) := (-l)"det 



ei 
DF^ 



(-l)"det 



ei 



Dip"" 



:= [J{Dip, Du), J^] xR, 

where DF^ G R" for f = 1, . . . n + 1, > is the jacobian of (p and (— 1)" 
is introduced to make sure the last component is positive. J {Dip, Du) is an 
M"-valued multilinear form, linear in the components Ui of Du and of weight 
n — 1 in the components of Dip. It is easy to check that J(I, Du) = —Du. 
The unit normal is: 

N{Dip, Du) = N{Dip, Du)/{\J{Dip,Du)\^ + (J^)^)^/^ 

Thus the angle condition may be stated in the form: 

/3[\J{Du,Dip)\^ + {J^^^l^ = J^|a^„, 

and we lose nothing by squaring it: 

B{Dip,Du) := PViDu,D^)\^ ^ PiW\9Do = 0- 



Du 



4. Compatibility and the choice of ipQ. Assume DipQ^Qj^^ = I. 
Differentiating in t the contact condition u^qdq = and evaluating at i = 0, 
we find: 

= c?'^(I, Duo)uoij = gl^UQij on ODq. 
To interpret this condition, consider the mean curvature at f = 0, on ODq: 

Ho = —[{J{I,Duo),gl^(poij) + J^od'^UQij], 

Vq 

where: 

vo = [\Jil,Duo)\' + = {\Duo\^ + = ^, 
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using (recall (3o := \/l — P^): 

J{I,Duo) = -Duq = -{DnUo)n = 
on ODq. Thus the compatibility condition is equivalent to: 

This implies we can't choose (po = Id (on all of Do), unless Hq^qu^ = 0, a 
constraint not present in the geometric problem (as seen above). Instead, 
regarding Hq as given (by Eq), and using: 

we find the compatibility constraint: 

ii^ij - PWn^)'Poij,n) = -^Ho on ODq. 

Po 

Given the zero and first order constraints on (pQ, this can also be written as: 

n'n^{(poij,n) = —^-^Hq on dD^. 
P Po 

The next lemma shows this can be solved. 

Lemma 4.1. Let Dq C be a uniformly C^+" domain (possibly 
unbounded), h e C°(5Do) (0 < o; < 1). 

(i) One may find a diffeomorphism ip G Dif P^°'{Dq) satisfying on ODq: 

ip = Id, d(p = 1, n ■ cP^p{n, n) = h. 

(ii) More generally, given a non- vanishing vector field e G C^^"(9-Doj K"), 
one may find tp G Dif f^~^°'{Do) satisfying on dDo: 

p> = Id, dnip = e, n ■ d^(p{n, n) = h. 

If BDq has two components, we may even require (p to satisfy the conditions 
in parts (i) and (ii) at diD^, O2-D0 (resp.), with different functions h. (This 
will be needed in section 11). 

As usual, a domain is 'uniformly (7^+"' if at each boundary point there 
are local charts to the upper half-space (of class C^+"^), defined on balls of 
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uniform radius, and with uniform bounds on the C^"*"" norms of the charts 
and their inverses. 

Remarks: 4.1. The proof is given in Appendix 1. 

4.2. Note that, in particular, ip satisfies the orthogonality conditions at 
dDo. 

4.3. It is at this step in the proof that we have a drop in regularity: 
for C'^'^" local solutions, we require C^^" initial data. While this is not 
unexpected in free-boundary problems (see e.g. P), I don't know a coun- 
terexample to the lemma if Dq is assumed to be a C^^° domain. 

4.4. In our application of the lemma, we in fact have h G C^~^'^{dDQ), 
but this does not imply higher regularity for ip. 

5. Linearization. The evolution equation and boundary conditions in 
'split gauge' are: 

' Ft-g'^{DF)Fij = 0, 
U\dDo = 0> 

B{D^,Du)\eDo = 0, 
[ 0{D^)^9no = 0, 

where: 

0{Dip) := {D^^,Dnip). 

Here D'^ip = Dip — {Dnip){-,n) is an R"-valued 1-form on ODq. We'll prove 
short-time existence for this system (with initial data 

by the usual fixed-point argument based on linear parabolic theory. Given 
F = [ip,u] in a suitable ball in this Holder space with center Fq = [ipQ,uo], 
it suffices to consider the 'pseudolinearization' of the system: 

Ft - (f \DFo)u,, = [(f^DF) - g'^{DFo)]F,, := .^(F,Fo); 

a fixed point of the map F ^ F corresponds to a solution of the quasilinear 
equation. 

For the nonlinear boundary conditions, we need the honest lineariza- 
tion at Fq. For the angle condition, a computation using the boundary 
constraints on uq and pQ yields: 

]^CoB[Dp, Du] = -f5f5oDnU - Pi{DnP, n). 

The corresponding linear boundary condition will be: 

(3(3oDnU + {I - (3^){DnP,n) = B{DF,DFq), 
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where: 

2B{DF, DFo) := B{D^, Du) - B{Duo, Dipo) - CoB[D{ip - ipo),D{u - uo)], 
and we used: 

~Co[D(po,Duo]\aDo = PPoDnUo + (1 - P'^){Dnipo,n)i9Do = 0. 

Also, B(DipQ, Duo)^a£)^^ = 0, so at a fixed point B{Dip,Du)\Q£i^^ = 0. 

Linearizing tlic orthogonality boundary condition, we find that CQO[Dip\ 
is the 1-form on SDq- 

(with sum over repeated indices.) The corresponding linear boundary con- 
dition is: 

{Dn^.prof) + {D^ip,n) = -n{D(p, Dipo), 

where: 

n{D<f, Dipo) := 0{D<p) - 0{Dipo) - CoO[Dip - Dip^], 
and we used: 

CoO[D^o\\dDo = {(DnV'of, •) + {D^m,n)\aD, = 0. 



6. Complementarity. We wish to apply linear existence theory to the 
system: 

Ft-g''iDFo)Fij=:F, 
with boundary conditions at dDo: 

u = 0, 

PPoDnU + Pi{DnV,n) =B, 

and initial conditions: 

Ut=0 = Uq, ipt=0 = V>0- 

It is easy to see that the initial data satisfy the linearized boundary condi- 
tions, and above we constructed ipo so as to guarantee g^-' {Duq, DpQ)uQij\Qj)^ = 
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0. (There is no first-order compatibility condition for ipo.) Thus the hnear 
system satisfies the required compatibility at i = 0. 

Since the linearized boundary conditions arc slightly non-standard, we 
must verify they satisfy the 'complementarity' (Lopatinski-Shapiro) con- 
ditions. We fix xq G ODq and introduce adapted coordinates (p, a) in a 
neighborhood A/o C A/' of in Dq: 



X 



ro + pn(cr), a={aa)eU, 



where Fq : ^ M" is a local chart for ODq at xq {U C W^~^ open). This 
defines a basis of tangential vector fields in ro(Z^), and we may assume 
that, at Xq: {Ta,Th) = 6ab and VTaTb{xo) = 0. Let U and tp be defined in 
(-pi,0) xWx [0,r] by: 

U {p, a, t) = u{ro{a) + pn{a), t), ip{p, a, t) = </?(ro(o-) + pn{a),t). 

In these coordinates, the induced metric is written (in 'block form'): 



[9] 



M + (Up? 



(V'p, V'a) + UpUa {tpa, Ipb) + UaU^ 







1 

In-1 



at t = and xq. 
We have: 



Upp = D^u{n,n) (since Vnn = 0), 



Uab = D'^u{Ta, n) + Du ■ V r^U = D^u{Ta, u) at Xq, 

and we don't need Upa, since gpa = at xq. 
Thus: 

trg,D^u{xQ) = p^D\{n,n)+J2D^<ra,Ta) = p'^Upp+J^Uaa := f3^Upp+A^U, 

a a 

and, likewise: 

trgoD^ip{xo) = li^Tppp + A^ip. 
For the linearized orthogonality operator, note that, at xq: 

Putting everything together, the linear system to consider at xq is: 



Ut - P^Upp - A^U = 0, 
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with boundary conditions: U\p=o = 0, 

Po{i'p,n) + (3Up\p=o = b{a,t), 
ii'p^Ta) + {tpa,n)\p=o = iOa{(J,t), a = 1, . . . n - 1. 

Laplace transform in t to 



Now take Fourier transform in cr € 
obtain: 



nn— 1 



,p G C,p < 0. 



In transformed variables, we obtain the system of linear ODE (in p < 0, for 
fixed {i,p)): 

(3^Upp-{p+\i?)U = 

(iHpp-{p+\i?)^ = Q- 
Writing the solution in the form: 



" U{p) ' 




- u{o) - 


_ m . 







we find the characteristic equation /3^7^ + p+ = 0, and choose the root 
7 so that ^7 = (1//3)\/A (where A = p -\- and we take the branch of ^: 

i2e(\/A) > 0). Here (p,^) G A, where: 

^ = {(p, G C X M"-i; \p\ + > 0, Re{p) > 

Thus the solutions decay as p — > — oo. Let be the space of such decaying 
solutions, dimcW"^ = n — 1. The relevant boundary operator on is: 



U 



U 

l3o{4>p,n)+(3Up 



\p=0 



11(0) 



(a vector in C x C x C""^). 

The 'complementarity condition' (see e.g. [3j) is the statement that 
B is a linear isomorphism from to C"^"*^. With respect to the basis 
{U{0), (Vi(0), n), {ip{0)Ta)} of W+, the matrix of B is (in 'block form'): 



[0]lx(n-l) 
[0]lx(n-l) 

[0](n-l)xl Ka](n-l)xl 



A 
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This is triangular with non-zero diagonal entries for every (p, ^) G A. Hence 
B is an isomorphism. 

7. Estimates in Holder spaces. 

For the fixed-point argument based on the linear system, we need esti- 
mates for ||.?^||a!, ||H||i+Q,, ||r2||i+Q,, of two types: 'mapping' and 'contraction' 
estimates. 

A bit more precisely, for T > 0, i? > and = Dq x [0, T] consider 
the open ball: 

(Fq = [(^0) '"o] is defined from the initial surface Sq, via Lemma 4.1.) Solving 
the linear system with 'right-hand side' defined by -F G defines a map 
F : F 1-^ F, and we need to verify that, for suitable choices of T and R, F 
maps into and is a contraction. 

Remark: The argument that follows is standard, and the experienced 
reader may want to skip to the statement of local existence at the end 
of the next section. On the other hand the result is not covered by any 
general theorem proved in detail in a reference known to the author, and 
some readers may find it useful to have all the details included. Another 
reason is that, although the 'right hand sides' are clearly quadratic, without 
explicit expressions one might run into trouble with compositions (which 
behave poorly in Holder spaces), or when appealing to 'Taylor remainder 
arguments' if the domain is not convex. 

For 'mapping', we need estimates of the form: 
1 1 J^(F, Fo) I la + I \B{DF, F»Fo) I + | |0(F>(^, D^o)\\i+a decays as T ^ 0+ , 
and for 'contraction': 

where //(T) ^ as T ^ 0+. 

Notation: The (a, a/2) norms are taken on Q-^, the (1 -|- a, (1 -|- a) 12) 
norms on ODq x [0, T]). Double bars without an index refer to the (2 + a, 1 + 
a/2) norm, single bars to supremum norms over Q^, and parabolic norms 
are indexed by their spatial regularity (a for {a, a/2), etc.) In general, we 
use brackets for Holder-type difference quotients. 
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We deal with the estimates for the 'forcing term' JT first. Consider the 
map 

g : Imm{W\W+^) GL^ 

which associates to the hnear immersion A the inverse matrix of {{Ai, Aj))^^-^^, 
inner products of the rows of A. Q is smooth, in particular locally Lipschitz 
in the space W of linear immersions. Hence if F^,F^ are maps — ^ R""^-*^, 
such that DF' G C°^'°/2(Q'r) and DF'{z) G K for all z G Q^, where K CW 
is a fixed compact set, we have the bound: 

\\g{DF') - giDF^)\\a < CKWDiF' - F^)\\^. 

In fact our maps F^ are in (72+a,i+Q!/2^ jjpi ^ p],Qj^ ^j^jg higher 

regularity we obtain the decay as T — ^ 0+. Assuming F^\t=o = -F^|t=o, we 
have: 

\D{F^ - F^)\ < [D{F^ - F2)]|^^rH^. 

To continue, we recall an elementary fact for Holder spaces: 

Let D C be a uniformly domain (not necessarily convex or 
bounded). Then if / G 0^(0) and a G (0, 1), we have: 

<Cz,||/|bi. 

Here 'uniformly C"*^' means D can be covered by countably many balls 
of a fixed radius, which arc domains of manifold-with-boundary local 
charts for D, with uniform bounds for the charts and their inverses. The 
constant Cd depends on those bounds. Applying the lemma to DF, where 
F = F^ — F^ vanishes identically at t = 0, and assuming T < 1: 

[DF]'^^ < c{\DF\ + \D'^F\) < c{[DF][^h^ + [D'^Fi^^T''/'^) < c||F||^"/^ 
(where c depends on Do) and similarly for the oscillation in t: 
[DFp^ < [L>F]|^VV2 < ||F||rV2, 

so we have: 

\\D{F^ - F'^)\\a < c\\F^ - F2||r°^/2. 
We conclude, under the assumption F^ = F^ at t = 0: 

\\g{DF^) - g{DF^)\\„ < ck\\F^ - F2||T"/2. 
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In particular, applying this to F and Fq, we find: 

WiGiDF) - g{DFo))D~^F\\^ < ck\\F - Fo||r"/2||F||, 
and for F^ and F^ coinciding at t = 0: 

\\{g{DF^) - g{DF^))D^F^\\a < ck\\F^ - F^WT'^/^WF^W, 
as well as: 

\mDF^) - g{DFo)){D''F^ - < ck\\F^ - Fo||r"/'||i^' - i^'ll, 

so we have the mapping and contraction estimates for !F{F, Fq) and J^{F^,Fo) — 

Lemma 7.1. Assume F,Fo,F^,F^ are in (72+"'1+°/2(qT. j^n+i^ 
have the same initial values, and that DF, DFq, DF^ , DF^ all take values 
in the compact subset K of /mm(M",R""'"^). Then: 

mF,Fo)\\a<CK\\F-Fo\\\\F\\T'^/^ 

mF\Fo) - J^iF'',Fo)\\^ < ck{\\F^\\ + - Fo||)r°/2||Fi - F^ 
In particular, if F G -B^: 

mF,Fo)\\a<coRT''/^ 

If F^,F'^ e B'^, we have: 

I \HF' , Fo) - HF^ , Fo) I U < coT°/2 \\F' - F\ 

(The constant cq depends only on the data at t = 0, and we assume T < 1, 
R < 1). 

Turning to the orthogonality boundary condition, first observe that: 

which has quadratic structure. Using a local frame (TaYaZl for TODq, we 
find the components O^: 

na{D^\ D</) = [A(^' - 'P^)Dj{^^ - ^o) + Dj{^^ - ip^)Di{ip^ - v?o)]nVi 
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(summation convention, i,j = l,... n), so Ha is a sum of terms of the form: 

b{x)D{ip^ — ip'^)D{(f)^ — Lp'^), where h{x) = n^T^ and the (p^ coincide at t = 0. 
It is then not hard to show that: 

\\h{x)D{^^ - cp')D{<p' - <^^)||i+« < c\\b\U+a\\^' - <p'\\\\cp^ - y^^IlT", 

with c depending on the norms of local charts for Dq. To bound the norm 
||n(8)ra||i+a, note \n\\Ta\ < 1, |D(n(g)Ta)| < |Dn| + |Dro| and [D{n'^Ta)]^^ < 
[Dnfx'^ + [Dto]". Since n = —{P/Po)Duo on dDo (and dDg is a level set of 
uo), we clearly have: 

\\Dn\\a + \\DTa\\a < c\\D'^Uo\\a < c\\uo\\. 

We summarize the conclusion in the following lemma. 

Lemma 7.2. Assume (f,(po G C^+"'^+"/^((5'^; R") have the same initial 
values. Then: 

||n(L»(^,L>(^o)||i+a < co||no||||^-(^o||'r" 

and 

\\n{D<p\D^^)\\i+a < colKIKII^' - ^o\\ + 11^' - ^o||)T"||^^ - ip\ 

with Co depending only on the data at t = 0. In particular, if F = [ip, u] G 
B^, we have: 

\\n{D<f,D^o)\\i+a<coR^T^, 
and for = [<f^,u^] e B]^, I = 1,2: 



To explain the estimates for the angle condition, we write the normal 
vector as a multilinear form on DF^: 

n+l 

N{DF) = Jn{DF) := (-1)" J](-l)^-^(i:>F^A. . . £>>^A. . . DF"+^)ej e 

i=l 

[DF'' omitted in the i*'*' term of the sum), where Z?F* € M" for i = 1, . . . , n+ 
1 and we identify the n- vector in R" with a scalar, using the standard volume 
form. The angle condition has the form: 

/?2|iV|2-(iV,e„+i)2 = 0on dDo, 
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and we set: 

B{DF) := (3^\Jn{DF)\'' - {Jn{DF),en+if, 
with linearization at DFq = [In\Duo\: 

CoB[DF] = 2p^{Jn{DFo),DJn{DFo)[DF])-2{Jn{DFo),en+i){DJn{DFo)[DF],en+i). 

Under the assumption F^ = F^ at t = 0, we need an estimate in 
^i+a,l±^ for: 

B{DF^, DF^) := B{DF^) - B{DF'^) - CqB[DF^ - DF^] 

= pWJn{DF^)\^ - \Jn{DF^)\^ - 2{Jn{DFo),DJn{DFo)[DF^ - DF^])) 

-{{Jn{DF^),el^^)^-{Jn{DF^),en+i?-2{Jn{DFo),en+i){DJn{DFo)[DF^-DF^^^ 

It win suffice to estimate the expression in the first parenthesis; the second 
is analogous. 

We need the following algebraic observation: if To = [Inl-Dtto] and T are 
n X (n + 1) matrices, the expression: 

|j„(ro + r)|2 - |j„(To)|2 - 2(j„(ro),Dj„(To)[r]) 

is a linear combination (with constant coefficients) of terms of the form: 

where the p(^2) (T) are polynomials in the entries of T (with constant coeffi- 
cients), with terms of degree: 2 < deg < 2n. 

Thus B{DF^,DF^) is a linear combination (with constant coefficients) 
of terms: 

uoip^2){DF^-DF^), uoiUojP^2){DF^ - DF^), p^2)iDF^ - DF^), 

with the p^2) as described; and hence is a linear combination of terms of the 
form: 

MF'J-F'jr, u,,u,,{Fl^-Fl^f, {Fl'-Fl^f 

(where 2 < d< 2n, 1 < j < n + 1,1 < i,l,k < n), which we write symboli- 
cally as: 

B{DF\DF^)^ b{x){DF^ - DF'^f, 

2<d<2n 
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where b{x) is constant or uoi{x) or uoi{x)uoj{x). For the degree d terms 
~ b{x){DF^ - DF'^Y, it is not hard to show the bound: 

||G'('^)||i+a < c||6||i+a||F' - fYt"", 2<d<2n. 
We conclude: 

Lemma 7.3.Assume F, Fq, F\ are in C^+°'''^+'^/'^{Q'^; M"+1) and have 
the same initial values. Then: 

\\B{DF,DFo)\\i+a < c(l + \\uo\\^)il + \\F- Fo||2«-2)r"||F - Fo\\\ 

\\B{DF\DF^)\\i+a < c(l + ||^^o||^)(l + - f2||2"-2)t"||F^ - F^, 
with c depending only on Fq. In particular, if F G B^: 

\\B{DF,DFo)\\i+a<coR''T'^, 

and if F\F^ e Sg: 

\\B{DF\DF^)\\,+a < coT^llF^ - F\ 
with Co depending only on Fq. 
8. Local existence. 

Let Dq C be a uniformly C^^" domain, not necessarily bounded 
or connected (note: we define our norms as the sum of the norms on each 
connected component). 

Given a C^"*"" graph Sq over Dq satisfying the contact and angle condi- 
tions, let tpo G Difp^°^ be a diffeomorphism given by lemma 4.1 (with the 
1-jet of the identity at ODq and 2-jet determined by the mean curvature of 
Eq at dDo). Then find uq G C^"'""(L'o) so that Fq = [ipQ^uo] parametrizes 
So over Dq. 

(Precisely, if [z, ^0(2;)] parametrizes So as a graph, and (^q is given by 
lemma 4.1, let uq = uqo ^q; so uo G C^"^"^.) 

We obtained in section 7 all the estimates needed for a fixed-point argu- 
ment in the set: 

B'^ = {Fe C^+'^'i+^/^^Q^^^n+i). - Foil < R,F\t=o = Fq}. 

Choose R < 1 and Tq < 1 small enough (depending only on Fq) so that, 
for F e B^j', F{t) = [ip{t),u{t)] defines an embedding of Dq, with ip{t) a 
diffeomorphism onto its image D{t). Let K C Imm{W,W+^) be a compact 
set containing DF{z) for all F G -Br, 2 G . Now consider T < Tq. 
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Given F G B^, solve the linear system (with initial data Fq) to ob- 
tain F G C^^"'^^"/^((5^)- (This is possible since the complementarity and 
compatibility conditions hold for the linear system.) This defines a map 
¥ -.F^ F. 

From linear parabolic theory (e.g. [3], thm VI.21]): 

\\F - Foil < M(||.F(F,Fo)|U + \\B{DF,DF^)\\i+^ + \\^{D^,D^^)\\i+^), 

where M > depends on the (7"'"/^ norm of the coefficients of the linear 
system, that is, ultimately on ||-Fo||. 

Prom lemmas 7.2-7.4 in section 7, it follows that: 

||F - Foil < Mco(i?T"/2 + i?2^") < R 

provided T is chosen small enough (depending only on Fo.) Thus F maps 
5^ to itself. 

Similarly, if F(F*) = F* for i = 1,2, standard estimates for the linear 
system solved by F^ — F^ give: 

||f1-F2|| <M(||.F(F\F2)||^^||^p^i^^^2)||^^^^pp-i^^-2)||^^^) 
Again the estimates in lemmas 7.2-7.4 imply: 

_ i?2|| < Mco(T"/2 + r")||Fi - F^ii < iiiJ?! _ i?2||^ 

assuming T is small enough (depending only on Fq). This concludes the 
argument for local existence. 

Theorem 8.1. Let Sq C M"+i be a (7^+" graph over F»o C W satisfy- 
ing the contact and angle conditions at ODq (Sq may be unbounded or not 
connected). There exists a parametrization Fq = [(/9o,^io] G C'^^°'{Dq) of Sq, 
r > depending only on Fo and a unique solution F G C^+"'-'^+"/^((5-^; M"^-'^) 
of the system: 

r dtF-g'i{DF)didjF = Q, 

with initial data Fo. For each t G [0, T), F{t) is a C^"*"" embedding parametriz- 
ing a surface Sj which satisfies the contact and angle conditions and moves 
by mean curvature. In addition, F{t) satisfies the orthogonality condition 
at dDf). 
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The hypersurfaces are graphs. For each t G [0,T), ip{t) : Dq — >■ D{t) 
is a diffeomorphism and T,t = graph{w{t)), for w{t) : D{t) ^ M given by 
w{t) = u{t) o ip-^{t). (We have w{t) G C'^+"\D{t)), 'less regular' than u{t) 
or <^(t).) -D(t) is a uniformly C^+" domain. 

Remark 8.1. This theorem does not address geometric uniqueness of the 
motion, given Sq. It only asserts uniqueness for solutions of the parametrized 
flow (including the orthogonality boundary condition) in the given regularity 
class. 

9. Rotational symmetry. In this section we record the equations for 
two rotationally symmetric instances of the problem: (i) Dq and D{t) are 
disks, and n > ('lens' case); (ii) -Dq and D{t) are complements of disks in 
M" ('exterior' case). For simplicity we restrict to n = 2. 

Let F(r) = [ip{r),u{r)] parametrize a hypersurface E, where ip{r) = 
(j){r)er is a diffeomorphism onto its image. Here 6^,651 arc orthonormal 
vectors, outward normal (rcsp. counterclockwise tangent) to the circles 
r=const. The unit upward normal vector and mean curvature are: 

1 



where: 



H = j^^^^^^{(l)rM{<l)r,Ur)[D\] - (u,e, , >i (<^, , Ur [I? V]) ) , 

M{(pr,Ur)[D'^u] = Urr + {<Pr +'U^)^^^^^, 
Mi<J)r,Ur)[D^^] = [At + (0? +«?)(^ - ^)]e^. 



Simplifying: 

1 u 
^ " (02+^2)3/2 i'l'rUrr - «r0rr + i4>r + ^r)^]- 

Now consider the time-dependent case F{r,t) = [4>{r,t)er,u{r,t)]. From 
the above expressions, one finds easily that the equation {dtF, N) = H takes 
the form: 

~ -lJ^r^-^i^r,y'r)[D'^u]) = Ur{er,<Pt - -T-^-^M{(j)r,Ur)[D^(p]) . 

(Py, ~r Uy, (py, Uy> 
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In 'split gauge', we consider the system: 



ut - ^ M{(t>r,Ur)[D'^u] = 0, 



Note that 0(r, t) = r solves the (f) equation, and that in this case the u 
equation becomes: 

wt - ^ = 0. 

l + w^ r 

This can be compared with the equation for curve networks: 

Wt - ^ = 0. 



The boundary conditions are easily stated (we assume Dq is the unit 
disk or its complement). 

The 'contact condition' at r = 1 is u = 0. For the 'angle condition' at 
r = 1, we find: 



Pi 



Assuming 0^ > at r = 1, this resolves as: 

Pur + Po4'r = at r = 1 (lens case); 

Pur — Po4>r = at r = 1 (exterior case). 

(For lenses, one also has at r = 0: u,. = and cpr = !•) Thus in both cases 
one can work with linear Dirichlet/Ncumann-typc boundary conditions. 

One reason to consider the exterior case is that (unlike the lens case) 
it admits stationary solutions. Geometrically, one just has to consider one- 
half of a catenoid, truncated at an appropriate height. For example, for 120 
degree junctions the equation for stationary solutions: 

H = m |r > 1|, 



1 + r 

Ur\r=l = VS, U\r=l = 



admits the explicit solution: 

V3, 



u{r) = -^(ln(2r + V4r2 - 3) - In 3), r > V3/2. 
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Problem. It would be interesting to consider the nonlinear dynamical 
stability of this solution (even linear stability is yet to be considered.) One 
may even work with bounded domains, by introducing a fixed boundary at 
some R> 1, intersecting the surface orthogonally (see Section 10.) 

10. Fixed supporting hypersurfaces. Extending the local existence 
theorem to the case of hypersurfaces intersecting a fixed hypersurface S 
orthogonally presents no essential difficulty. The case of vertical support 
surface leads directly to graph evolution with a standard Neumann condition 

on a fixed boundary; we consider the complementary case where 5 is a graph. 
Let S C M""''"'^ be a embedded hypersurface (not necessarily connected), 
the graph over V C M" of S G C^(P), oriented by the upward unit normal: 

Ky) :=—%), Hy)-=[-DB{y),l]eR^xR, vb := V^T\DB(^ . 

Vb 

V is assumed to be nowhere vertical in T) (DB ^ 0). To state the problem 
in the graph paramctrization, we consider a time-dependent domain D{t) C 
MP with boundary consisting of two components d\D{t) and d-zDit), both 
moving. The hypersurface St is the graph of w{-,t) over D{t), solving the 
parabolic equation: 

wt - g'\Dw)wij = in := (J D{t) x {t} G x [0, T], 

te[o,T] 

with boundary conditions: 

0|9iD(t) = 0, v'l + I^W^P|Sii3(t) = V/5 

(as before), and on dzDit): 

w = B, Vw VB = -1. 

(The first-order condition on d2D{t) is equivalent to (i^, A^) = 0). 

Differentiating in t the boundary condition w = B leads easily to an 
equation for the normal velocity of the interface V{t) = d2D{t): 

. vH 

Bn - Wn 

Note that Wn at d2D{t) can be computed from Bn, since: 
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in particular neither i?„ nor Wn can vanish (so both have constant sign 
on connected components of 52-D), and one easily computes: Wn — Bn = 

-Vl/Bn. 

Let A = S n 5 be the {n — l)-manifold of intersection, the graph of w 
(or B) over d2D. Given the graph parametrizations of S and S: 

G{y) = [yMy)\, B(y) = [y,B(y)], y G 
and T G Td2D, we have the tangent vectors: 

Gn ■■= [n, Wn] G TS, Gb ■■= [VB, -1] = -vsiy G TE, Gr := [r, V-uj-r] G TA, 
and the second fundamental forms of S and 5 (for e G arbitrary): 

A{dGe, dGe) = -d'^w{e, e), ^(dBe, cffie) = — <i^S(e, e). 

V VB 

From {i^,N) = at d2D, it follows easily that (cp. [9j): 

A{Gr, v) = -A{Gr, N), T G TdD. 

For the remainder of this section, we concentrate on the boundary con- 
ditions at (92 -Do, and denote this boundary component simply by ODq. To 
establish short-time existence, we consider as before the parametrized flow: 

Ft-trgd^F = 0, g = g{dF) , F=[ip,u]. 

The contact and angle boundary conditions are: 

Again we have two scalar boundary conditions for n + 1 components. Here 
the solution is easier than at the junction. With the notation Fn = dFn = 
[iPn,Un], we replace the angle condition by the 'vector Neumann condition': 

Fn _L TS, or Fn = —avsv on SDq, 

where a : dD^ — > M, or equivalently (since this leads to q = —Un)- 

= -Un{VB o Lf) on (92-Do. 
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Clearly the Neumann condition implies the angle condition {N, v o ip) = 

0, but not conversely. This linear Neumann-type condition can easily be 
incorporated into the fixed-point existence scheme described earlier. 

There is one issue to consider: the and Ist-ordcr compatibility con- 
ditions must hold at ODq, at t = 0. The initial hypersurface Sq uniquely 
determines wq and Dq C (satisfying wq = B and Vwq ■ VB = — 1 on 
dDo), and then once tpo G Diff{Do) is fixed, uq = woocpo is also determined. 
We may assume: 

(fio = id, (fon = V-B on dDo, 

so: 

Uon = Vwo ■ fOn = VlfQ • VB = - 1 on ODq, 

and then the Neumann condition i^on|92£'o ~ ~'^^-B^ holds at t = 0, on ODq. 
The first-order compatibility condition is: 

tVgd^UQ = ut = VB ■ (ft = VB ■ tVgSipQ on dD^, 

or equivalently: 

trg{iy,d^Fo) = on ODq. 

(This is not a mean curvature condition; the mean curvature of So is = 
trg{N,cPFo).) 

From now on we omit the subscript 0, but continue to discuss compatibil- 
ity at t = 0. First observe that the Neumann condition leads to a splitting of 
the induced metric. Given r G TODq, let F-^ = dFr G TA. Then (recalling 
Un = —1 on dDo): 

{Fr, Fn) = {[t, dBr], [ipn, Un]) = V B ■ T ~ V B ■ T = 0. 

Thus we have: 

trg{iy,d'F) = g-\i,,d^FiTa,n)) + g''''W,d^F{Fn,Fn)), 
for a local basis {T^ = dFTaj^Zl of TA, with gab = {Ta, T^) and 5„„ = 

\Fn? = vl. 

Differentiating in n the condition Un = Vw ■ (pn ( assuming, as usual, n 
extended to a tubular neighborhood M of dDo as a self-parallel vector field) , 
we find: 

Unn = d^w{n, VB) + Vw ■ d^(p{n, n). 
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This is used to compute: 

{v,cPF{n,n)) = — K„ - VB ■ (f(p{n,n)] 

VB 

= —[(fw{n, VB) + {Vw - VB) ■ (f(p{n, n)] 

VB 

= -vA{Gn, 1^) + —{wn - Bn)n ■ cf(f{n, n). 

VB 

Bearing in mind the expression for Wn — B^ found earher, the compatibihty 
condition may be stated in the form: 

^n-d\{n,n) = -vA{Gn,u) + g-\SF{Ta,n),v). 

We are now in the same situation as in section 4: given the 1-jet of ipo on 
ODq , we extend ipo to a tubular neighborhood M of SDq (and then to all of 
Do)) so that n-d?ip{n, n) has on ODq the value dictated by the compatibility 
condition (using Lemma 4.1(ii)). We just need to verify that the right-hand 
side of the above expression depends only on Sqi 5, and the 1-jet of (po over 
dDo. Clearly only the term g°''^{i',(fF{Ta,Tb)) is potentially an issue. 

Fix p G ODq, and let {tq} be an orthonormal frame for TODq near p, 
parallel at p for the connection induced on ODq from R". If /C denotes the 
second fundamental form of BDq in M", we have: 

Ta{n) = f^{ra, n)n ( at p) 

(on the left-hand-side, r;, is regarded as a vector- valued function in M"). Still 
computing at p, this implies: 

(fF{ra,n) = TaidFn) - dF{ra{n)) 

= Ta{dMTb) - /C(Ta, rb)F„ 
= d^B(ra, n) + /C(Ta, n)Mn - }C{Ta, Tn)Fn, 

where F„ = —vu and B„ = dMn G TS. Hence: 

{U, d'FiTa, n)) = {U, d^M{Ta, Tfe)) + vJC{Ta, u) = A{Ta, Tfe) + vJC{Ta, u). 

This clearly depends only on S and on Sq. We summarize the discussion in 
a lemma. 
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Lemma 10.1 Let Sq = graph{wo) be a graph over Dq C (a 
uniformly domain), intersecting a fixed hypersurface S = graph{B) over 
ODq. Consider the parametrized mean curvature motion with Neumann 
boundary condition: 

F e C^^\Do X [0,T]) ^ M.''+\ F = [ip,u] 

Ft - trg(fF = 0, g = g{dF), uo ip = B and Fn -L TS on ODq. 

Then (po G Diff{Do) can be chosen so that (with uq = woo tp^) the initial 
data Fo = [(po , uo] satisfies the order zero and the first-order compatibility 
conditions at ^ = and ODq: 

Von = -non(VS o (pq), {u o ipQ, trggcfFo) = 0. 

Remark 10.1. Differentiating dwTa = dBra along Tf,, we find: 

d^w{Ta, n) - d^B{Ta, u) = {Wn - Bn)IC{Ta, u) 

(reminding us that, although w = B on ODq, the tangential components of 
their Hessians do not coincide.) From this follows the expression for IC in 
terms of A and A: 

Hra^n) = ^—^[vA{Ta,n) - VBAiTa,n)]. 

Wn — Bji 

It is also easy to express the corresponding traces in terms of the mean 
curvatures and of A in S and S: 

H'' = —g-'A{Ta,n), n'' = ^g''''A{Ta,n). 



11. A continuation criterion. Once local existence has been estab- 
lished, it is easier to obtain geometric estimates (in particular using the 
maximum principle) for the solution in the graph parametrization. (From 
this point on , we focus on the 'lens' case, without fixed support hypersur- 
faces.) 

For a time interval / = {to,ti) C [0,T] set: 

E = {z={y,t)eW'xI;ye D{t)}, S = {{y, ty,t e I,y e dD{t)}. 
Let w he a solution in E of: 

wt - g'^ {Dw)DIjW = 
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with boundary conditions on S: 

w = 0, DnW = Pov, ■= \/l - /?^- 

For the remainder of the paper we assume D{t) is bounded, for each t £ I. n 
denotes the (t-dependent) inner unit normal at dD{t), extended to a C^'^ 
unit vector field in a tubular neighborhood of D{t) so that Dnn = 0. 

Denote by L the operator L = dt — g^^ {Dw)didj , so Lw = in E. The 
following height bound is immediate. 

Lemma 11.1. Assume < wq < M in D^Iq). (If there is a support 
surface S, we assume B{y) > in D and M < sup-p B.) Then < w < M 
in E. 

Proof. Follows from the maximum principle applied to L, since < w < 
M holds on the parabolic boundary dpE. 

It is well-known that the function v = y^l + \Dw\^ solves the evolution 
equation (assuming Dw S C'^'^{E), see e.g. ^): 

Lv + '^g'^ViWj = -v\A\l. 

From the maximum principle, we have the following global bound on v 
(equivalently, on 

Lemma 11.2 Assume w is a solution with Dw G C'^'^{E). Then we 
have on ^ 

v{z) < max{sup£,(io)u(2;,to), 

Proof. By the maximum principle, max^v = maxg^^ju. Note v^g = ^. 

It follows from this lemma that gij{t) is uniformly equivalent to the 
euclidean metric in D{t): if < "u in E, and X is a vector field in D{t): 

\X\l < \X\l = g^JX'X^ = \X\l + {X ■ Dwf < \X\l{l + \Dw\^) < v^\X\l. 

Also, if oj := v~^Dw: 

yZ yZ yZ 

This equivalence of norms clearly extends to tensors, in particular to h: 

-\h\l<\h\l<Cn\h\l 
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where (throughout this section) Cn denotes a constant depending only on n 
and V. More generally, defining: 

ml := J2(9kh^J)^ ml ■.= ^ g'^'iOkh^jmh^j), 

i,j,k i,j,k.,l 

we have, for each pair z, j: 

E(^fc^ii)' = ^6'''dkhijdihij = ^{g'''+M)dkhijdihij 

k kjl k^l 



1 

= {d^hijf + ^g^^dkhijdihij < ^g^^dkhijdihij + (1 - 3 ) ^{dkhijf, 

k -.Jj k ^If k 

and hence, adding over i, j we have: 

mil < mil < v'mil 

The same argument works for second derivatives. The norms defined by: 

\d^h\l := ^ {dmdkhijf, \d^h\l := ^ g""\dmdkhi,){dndkhij) 

i,j,k,m i,j,k,m,n 
are uniformly equivalent in E: 



The point is that these euclidean norms satisfy easily computed evolution 
equations. Using the results in Appendix 2, it is straightforward to see that: 

L[\h\l] = -2\dh\l + 2Y,Cijhij, dj := L[hij], 



L[\dh\l] = -2\d'h\l+2j2{dkQj){dkh^,)+2 J2 idkg"'nidmdnhj){dkh,j) 

i,j,k i,j,k,m,n 

In symbolic notation, we have: 

Cij ^ h * h * h, dkCij (dh) * h* h, 
which combined with the previous remarks implies: 
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for a constant c„ as above. In addition, from (see Appendix 2): 

9^5"^" = hJ^oj'' + hloj"^ and \g'^\ < 2, < 1, 
we have ISfcS™'"! ^ Cn|/i|g for each m,n,k. We conclude: 

L[|/i|^]<-2|a/.|2 + c.|/.|^, 

L[\dh\l] < -2\d^h\l + cn\dh\l\h\l + Cn\h\g\d^h\g\dh\g. 

These differential inequalities imply the continuation criterion given in Propo- 
sition 11.3. 

Recall that for mean curvature flow (or mean curvature motion) of 
graphs, interior estimates for v imply interior estimates for A and its co- 
variant derivatives V™A of any order (see [1] or [2]). In the following 'con- 
tinuation criterion', global bounds are needed. 

Proposition 11.3. Assuming T^ax is finite, let w : E'^^°-'^ ^ M be a 
maximal solution, defined for t G [0,Tmaa;)- Then: 

limsup( sup -|- sup \V h\ g{y, t)) = oo. 
t^T (y,t)eE yedD{t) 

Proof. By contradiction, assume we have bounds in [to^^i] (for ti arbi- 
trarily close to Tmax)'- 

sup|/i|g<ao, sup \Vh\g<bo. 
zeE 2eSuD(o) 

For a > to be chosen (small), define the function on E: 

f{x,t) = a\dh\l + \h\l 

Then, for any > 0: 

L[f] < -2a\d'^h\l - 2\dh\l + Cni4 + aal\dh\l + aao\d'^h\g\dh\g) 

< -2a|9^/i|^ - 2\dh\'^g + Cnaoar]\d'^ h\'^g + {cnCtal + '^"^°" )|9/i-|g + CuOq. 

Choosing rj so that CnOor] < 1, then a so that c„aao -|- '^"^^ < 1, we ensure 
that: 

L[f- Cn4t] < 



31 



in E. By the maximum principle: 



a sup \dh\1 < sup / < sup / + CnagT < Cn{aQ + abg + a^T). 

E E dpE 

This imphes a uniform C^{D{t)) bound for w in E, and hence (by hnear 
paraboUc theory, given the uniform bound on \ Dw\ from lemma 10.2) a C^"*"" 
bound for some < a < 1. So we can apply the local existence theorem 
with initial data Sj^ , to continue the solution for a time depending only on 
bounds at to, contradicting the maximality of Tmax- 

Lemma 12.2 (in the next section) implies the conclusion can be strength- 
ened: only a uniform bound on tangential covariant derivatives of the second 
fundamental form /C of the moving boundary (in M") is needed: 

Proposition 11.4. Assuming T^ax is finite: 

limsup[ sup + sup \V rK.\g{y,t)\ = CO. 
t^Traax (y,t)(iE yedD{t) 

It is possible to strengthen this further and show that: 

limsup[ sup \h\g{y) + sup \VrlC\g{y,t)] = oo. 
t^Trnax yedD{t) yedD{t) 

That is, the supremum of \h\g on the moving boundary controls its value 
in the interior. The reason is that we already have a bound on sup^v; as 
remarked earlier, it is a well known-fact for mean curvature flow of graphs 
that this implies interior bounds for the second fundamental form and its 
covariant derivatives ([1], [2]). In the next lemma we describe a global 
argument for mean curvature motion of graphs with moving boundaries. 

Proposition 11.5. Let w : E ^ M he a solution of graph m.c.m in a 
spacetime domain E C x [0, T], where T < oo. Assume the first derivative 
bound v{x, t) < v holds globally in E. Then if the bound \h\g < /iq holds on 
the parabolic boundary dpE, we also have the global bound: 

\h\g < ao in E, 

for a constant oq depending only on n, v, ho, T and the initial data of w. 
Proof. The idea is to consider a function in E of the form: 

ip = AvP + \h\y + B\h\l, 
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where A, B and p are positive constants. We claim it is possible to choose 
these constants, and also C > (all depending only on n and v) so that: 

L[(p] <CmE. 

Thus L[ip — Ct] < in and hence by the maximum principle: 
supS|/i|g < sup V? < supo^E^P + CT, 

E E 

which clearly implies the bound claimed in the proposition. 

The proof that ip as above exists is (of course) based on the evolution 
equations for \h\g and v (see Appendix 2), which imply (for constants Cn, dn 
depending only on n): 

L[\h\l]<-2\Vh\l + Cn\h\% 

L[vP\ = -pvP\h\l - p{p - l)vP-^\dv\l - 2pvP-^g'^ViWj. 
Here \dv\g := g^^v^vi, and we have the bounds: 

L[vP] < -pvP\h\l-p{p-l)vP-^\dv\l + dnPvP-'\dv\g, 

L[vP] < -pv^lhff, - \p{p - 1) - \]\dv\l + dlp\vf^P-'\ 
The main term in L[ip] is: 

=L[|^|>f+LK]|/i|2-4/V-'5fe«(/i,V,/i)g := (7)+(77)+(m), 

where: 

{I)<-2vP\Vh\l + CnvP\h\l; 

{III) < 4pvP-'\dv\g\h\g\Vh\g < ^ | V/l | ^ + 7^ (^ - 1 ) 1;^" ^ | ^ | ^ , 

for an arbitrary constant 7 G (0, 1). With 77 > to be chosen sufHciently 
small later, we estimate the last term in (II): 

dnPvP-^\dv\g\h\l < vdnPvP-^\dv\l\h\l + ^dnPvP\h\l. 

Adding to these estimates for (I) + (II) + (III) the term L[B|/i|^], we have: 

L[\h\y + B\h\l] < [-^ - 2 - ^^BV\Vh\l + [c^-p+ '^;BV\h\t 
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+ [r]dnP - (1 - l)p{p - l)]vP ^ 



\dv\l\h\l + [^dnPV\h\l. 



Given 7 € (0, 1) arbitrary, we choose p > so large that (2/7)c„ < p — 1, 
then i? > so that: 



In this way we ensure that, in the expression above, the coefficients in the 
first two square brackets are negative. Choosing r] > sufficiently small 
(depending on 7 and p), the same holds for the third square bracket. Finally, 
in view of the second estimate given above for if we add ^[^41;^] with 

A > {2/rj)dn we also take care of the last square bracket (we also assume 
p{p — 1) > 1/4), and then: 



concluding the proof. 

12. Boundary conditions for the second fundamental form. 

Proving global existence for the mean curvature motion of graphs over 
time-dependent domains requires estimates for the second fundamental form. 

The simplest form of the evolution equations for hij and H = g^^hij is 
given in terms of the differential operator on functions: L[f] = dtf — trgd^ f . 

The evolution equations for hij and H are given in Appendix 2. In 
this section we derive boundary conditions for h and H; the development is 
similar to work of A. Stahl [9] for MCF of hypersurfaces intersecting a fixed 
boundary orthogonally. 

It is easy to see that h splits on dD{t): if r G TdD{t) is a tangential 
vector field, and n = is the inner unit normal: 

hin, t) = -d'^win, r) = -(r(tt;„) — Dw ■ Vru) = on dD(t), 

V V 

since Wn = Po/P on the boundary and Vi-n G TdD{t) (V is the euclidean 
connection.) In particular, it follows that h{Dw,T) = on dD{t). 

Boundary condition for H. We derived in section 2 the equation for the 
normal velocity of the moving boundary Tt = dD{t). Letting T{0,t), 9 G 
be any parametrization of Tt, we find for r„ := dtT ■ n: 



2p 




L[AvP + \h\y + B\h\l] <C := Adlp\vf^P-^^ 



tn = H = -—H Sit dD{t). 
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Since {N,en+i){t,T{t)) = /? on dD{t) we have: 

where d^N = —g^^hi^Gj, with e„+i component: 

{dkN,en+i) = -g'^Wjhik = — \h{Dw,dk) = — \wnh{n,dk). 

V V 

Hence we find, on dD{t): 

{dtN,en+i) = ^h{n,t) = ^tnh{n,n) = -fiHhnn- 

On the other hand, using dtN = — V^-ff — Hv'^V^v, combined with the 
expressions (vahd on dD{t)): 

(V^i^,e„+i) = g'^Hi{Gj,en+i) = g'^H^wj = \wiH, = = f3PoHn, 

2 

\V V,en-\-i) — 7^ — TTilnn — Pn'^nni 

yZ yZ 

we find on dD{t): 

{dtN, Sn+l) = -m{Hn + PoHhnn). 

Comparing these two expressions for {dtN,en+i) yields: 

m 

a Neumann- type condition for H on dD(t). 

Boundary conditions for hij. Fix p G dD{t) and let (tq) be an orthonor- 
mal frame for TpdD(t) (in the induced metric), satisfying V^^Tb{p) = (V'" 
is the connection induced on by V, or, equivalently, by V, the Levi- 
Civita connection of the metric g in D{t)); we extend the Ta to a tubular 
neighborhood so that V„ra = 0. Differentiating h{n, r^) = along r^, we 
find: 

{Vr,h){n,Tb) = -h{Vr,n,n) - h{n,Vr,n). 

The second fundamental form /C(r, r') of T in (D{t),eucl) (equivalently, in 
{D{t),g)) is defined by: 

VraTb = Vl^n + IC{Ta, n)n on dD{t). 
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To relate K, to h\dD(t)-, iiote that since = on dD{t): 

h{Ta,n) = {[Vr,n,o],N) = -Vr^n ■ ^ = -PolC{Ta,n). 

(So we see that Tt convex with respect to n corresponds to St concave over 
D{t), as expected). In the appendix we observe that Vg.dj = {hij /v)Dw. 
Then: 

= Tb + IC{Ta,Tb)n + —h{Ta,Tb)n = {- + Po)h{Ta, n)n = --^h{Ta,Tb)n 

V Po PO 

at p, given our assumption V^^Tb{p) = 0. We use this immediately to 
compute, at p: 

^2 I 

^Ta^ = {^Tan,Tb)gn = - ijl ,V raTb) gU = — | U | ^ /l (t^ , Tfe ) Tfe = —h{Ta,Tb)Tb, 

PO PO 

since |n|^ = gijU^n^ = 1 + = /?~^ at p. We conclude, using the Codazzi 
equations: 

1 

(V„/l)(r„,rfc) = {Vr^h){n,Tb) = -^y^h{Ta,Tc)h{Tc,Tb) + —h{Ta,Tb)hnn- 

PO ^ PO 

This can also be written in the form: 

/3o(V„/i)(r,r') = -{h'-^f{T,T')+P'Kr.HT,T'). 

It turns out the expression for covariant derivative of h with respect to the 
euclidean connection V is exactly the same (at dD{t)): 

/3o(V„/i)(r,T') = -(/i*"")2(r,T') +/32/i„„/i(T,r'). 
The reason is that VnTa = at the boundary, also for the g-connection: 

VnTa = VnTa + nWlV Q^dj = + ^h{n, Ta)Dw = 0, 

so in fact: 

{Vnh){Ta,Tb) = n{h{Ta,Tb)) = (V nh) {Ta , u) . 

As done in [9|, we combine this with the result for Hn to compute 
(Vnh){n,n). From: 

Hn = Vnitvgh) = trgiVnh) = /J^ (V„/i) (n, n) + ^ (V„/i) (t„ , , 
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we find: 



Po Po Po 

1 tan\2 I Mi,2 \ _ ^ 1 1,|2 



Po Po 



since 5"" = at dD{t). Equivalently: 

1 



/3o(V„/i)(n,n) = on aZ)(i). 



It is easy to obtain the corresponding expression for the euclidean con- 
nection. Noting that at dD{t): 

V„n = V„n + n^n^-hijDw = fS^hnnn, 

V 

we find: 

(ynh){n, n) = n{hnn) = {Vnh){n, n) + 2/i(V„n, n) = {S/nh){n, n) + 2/?o/i^„, 
so that: 

/3o(V„/i)(n,n) = + 2/32/iL on aL>(t). 



It turns out that the expressions just derived, combined with the max- 
imum principle proved in are not enough to estabhsh that concavity is 
preserved. We derive a suitable maximum principle in section 13. The result 
of the next lemma yields a continuation criterion stated earlier (Prop. 11.4). 



Lemma 12.2. Let w{y,t) be a solution of graph MCM, with constant- 
angle boundary conditions, in E G ffi" x [0, T). Denote by /C the second 
fundamental form of Tf = dD{t) in R". Suppose that, for some oq > 0: 

sup{|^|(y,t) + |V,/C|(y,t);y G dD{t),t e [0,T),t G TydD{t),\T\ = 1} < cq. 

Then also: 

sup{|VA|;y E dD{t),t e [0,T)} < 00. 

Proof. From the boundary conditions computed above for V/i, we have 
at boundary points: 

|(V„/i)(T,r)| + |(V„/i)(n,n)| + |(V,/i)(n, t)| < cq, 
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where cq depends only on (3 and oq. The remaining components of Vh are: 

{yTh){n,n) = {Wnh){T,n) and {VT'h){T,T), 

and since /i*"" = —(^qJC at boundary points, the last one is assumed bounded 
in [0,r]. In addition, on dD{t): 

t{H) = (3\Vrh){n, n) + 20'h{Vrn, n) + 5][(V^M(t„, + 2/i(V^r„, r„)], 

a 

with all terms on the right bounded, except for the first one. Thus a bound 
on (Vr^)(n, n) would follow from a bound on r(i/). But this follows from the 
uniform gradient estimates (up to the boundary) of linear parabolic theory, 
since is a solution of (see Appendix 2; uj := Dw/v): 

dtH - trgCpH = \h\lH + Hh^{u,u) - H^h{u,u), Hn\dD{t) = ^Hhnn, 

in which all the coefficients are uniformly bounded in [0, T). The bound 
depends only on gq and the initial data. (The hypotheses of the proposition 
imply that the necessary regularity conditions on dE are satisfied.) 

Finite existence time. 

In the next section we show that weak concavity at t = is preserved 
by the evolution. Assuming this, it is not difficult to derive that the flow is 
defined only for finite time. 

Lemma 12.4. Let w{y,t), {y,t) e E cWx [0,T) define a graph MCM 
Tit with constant-angle boundary conditions on a moving boundary. Assume 
So (and hence S^, for all t) is weakly concave. Then: 

Assume i?|t=o < Hq < (where Hq is a negative constant). Then 
T < t^ = ^rr^ ( we are assuming T = sup{t G [0,T);D(t) ^ 0}). Here 
Cn > depends only on n and an upper bound for v m E. 

The proof is based on the evolution equation and boundary condition 
for H (see Appendix 2: iv = Dw/v): 

L[H] = \h\lH + Hh?{u,u) - H^h{u,u), = {0" / (3^)Hhnn. 

Since h'^{oj,oj) > 0, \h\l > {l/n)H'^ and (given that /t < 0) h{oj,oj) > 
\Dw^H, we have: 

L[H\ < -H^ + \Dw\'^H^ < CnH^, 
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where Cn depends on n and on sup^; \v\ (already known to be finite). Let 
^(t) solve the o.d.e. = c„0^,^(O) = Hq: 

cp{t) = Ho[l - 2cnH^t]-^/^, 0<t<U:= — 

Then with := {l/n){H^ + H(f) + 4?) > 0, setting x = H - (f): 

L[x] < V'X in E; 

Xn = -rU + ^^Kn > -rX on diE 

Po PO 

(since ^ < and hnn < 0). Given that x < at t = 0, it follows from the 
maximum principle that x < 0, or iJ < ^ in [0, min{T, t*}). This shows 
< T is impossible, since cf) ^ —oo as t ^ t*. 

Remark 12.1. It would be natural to try to show that a negative upper 
bound Hq on the mean curvature (at t = 0) is preserved, at least under 
the assumption of concavity. Unfortunately, the evolution equation for H 
(under graph m.c.m.) does not lend itself to a maximum principle argument. 
Letting u := H — Hq, we have 



L[u] = \h\lu + uh^{uj,u)) - u{H + Ho)h{LO,uj) + HqQ in E, 

Q:=\h\l + h^{uj,uj)-Hoh{u,u). 

At a point where n = 0, we would need to show L[u\ < 0. But it is not 
true that Q > at such a point, even when n = 2. (?x„ > does hold at 
boundary points.) 

13. A maximum principle for symmetric 2-tensors. 

In this section we prove a weak maximum principle for the parabolic evo- 
lution of symmetric two-tensors on bounded euclidean domains, with moving 
boundaries and Neumann-type boundary conditions. The hypotheses are as 
follows. 

Let E C M^x [0, T] be connected, open and bounded (with boundary), 
with D{t) = En (M" X {t}) bounded, open, connected for each t G [0,r]. 
Fix > so that E C Br{0) x [0, T]. 

On the 'lateral boundary' of E: 

diE := {z = {x,t);t G [0,r],a; G dD{t)}, 
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we define the inner unit normal n = nt E M". Extend nt to a vector 
field in all of D{t) (so that it is in C^'^(£^, M"), arbitrarily except for the 
requirements that |n| < 1 pointwise and V„n = in a tubular neighborhood 
of dD(t). (Here V denotes the euclidean connection, so this requirement can 
be written n^diU^ = for each j.) 

The assumptions on the coefficients are given next. 

g = gt 'is & f-dependent Riemannian metric in D{t), uniformly equivalent 

to the euclidean metric for t G [0, T] ; 

X = Xt'is & bounded t-dependent vector field in D{t), satisfying X-ra > 
for z G diE; 

q = q{z,m) assigns to each z E E and each m in S (the space of quadratic 
forms in R") a quadratic form q £ S. q is assumed to be C^'^ in z, locally 
Lipschitz in m (uniformly in z £ E); 

b = b{z, m) G S is defined for z G diE, with the same regularity assump- 
tions. 

Theorem 13.1. Assume m G C^'^(£^;§) satisfies in E the differential 
inequality: 

dtuiij - trgCpruij < X ■ drriij ■ +qij{-, m{-)), 
and on diE the boundary condition: 

n ■ dmij{z) > bij{z, m{z)). 

Suppose the functions q and b satisfy the following 'null eigenvector condi- 
tions': if, for some m G §, 1^ G M" is a null eigenvector of m {rhijV^ = OVi), 
then, for any z e E (resp. any z G diE): 

qij{z,m)V'V^ < (resp. bij{z,m)V'V^ > 0). 

Then weak concavity of m at f = is preserved: 

m < in D{0) ^ m < in ^. 



Proof. The assumptions imply there is K > (depending only on E and 
on the functions X, g, n, q and b) satisfying: 

\n\c^AiE)<K, \X{z)\eucl < K, \g{z)\ + \g-\z)\<K, z e E; 
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and ii m,m e C'^'^{E,S) satisfy (for some fx : E ^ M+): 

—li{z)l < m{z) — m{z) < n{z)I 

(where I = (Sij) and the inequahty of quadratic forms has the usual mean- 
ing), then also: 

q{z, m{z)) < q{z, m{z)) + KiJ,(z)I, z & E, 
b{z, m{z)) > b{z, m{z)) - KiJ.{z)I, z e diE. 

Now define, ioi z e E: 

ip{z) := -2Kn{z) ■ x := 2Ks{z), 

where we use the euclidcan inner product and, on diE, s is the 'support 
function' of dD{t) (positive if D{t) is convex and contains the origin). It is 
clear we may find M = M{R,K) > depending only on K,R and |n|(^2,i 
so that: 

\ip\c2,i < M, \dtp\l + \trgd^ip\ <M, \X ■ dip\ < M. 

We assume also M > K. Now, given m as in the statement of the theorem 
and given constants e > 0, 7 > and 6 > 0, define for z E E^ := Er]{t < 6}: 

m{z) := m{z) - (et + 7e^(^))I, z G E^. 

Clearly m G C^^^(E^;§). We now derive the constraints on 6, e and 7. It 
will turn out that 6 must be taken small enough (depending only on K, R), 
e > is arbitrary and 7 is e times a constant depending only on K, R. 

The following inequalities arc easily derived: 

q{z, m{z)) < q{z, 7n{z)) + K{et + 7e'^(^))I; 

X-dm = X-dm + j{e'^X ■ d(p)I <X-dm+ {je'^M)I; 
dtm = dtm - el - {^e'^dtif)! < dtm + {je'^M)I - el; 
trgd^rh = trgd^m - 76*^(1^(^1^ + trgd^ip)! > trgd^m - {-fe'^M)!. 

We use this to compute: 

dtm - trgd^m < dtm - tVgd^m + (27e^M)I - el 
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< q{z, m{z)) +X-dm + {2-fe'^M)I - el 
< q{z, m{z)) +X-drh + K{et + 7e'^)I + (3M7e'^)I 
< q{z, m{z)) +X-drh + Metl + AMjeH - el. 
We conclude the inequahty: 

dtifi — tTgd^rh < q{z, m{z)) + X ■ dm — (e/2)I 

will hold in E^, provided the constants are selected so that, for z e E^: 

4M7e'^(^) + Met < e/2. {A) 

Turning to boundary points z = (x, t) G diE, note that dn^ = — 2-fC, so 
that: 

dnm{z) = dnm{z) - (7e^(^)d„^(^))I > h{z,m{z)) - {^e"^^'^ dn^{z))l 

> b{z, m(z)) - K{et + 7e'^(^))I - {je'^^^Un(p{z))I 
> b{z,m{z)) + K(7e'^W - et)I, 

so that the inequality: 

dnrh{z) > b{z,m), z G diE^ , 

will hold provided the constants are chosen so that, on diE^: 

et<je'^^'\ (B) 

Bearing in mind that, on E: e^^^^ < e'^^^'^ < e^^^, it is not hard to arrange 
for (A) and (B) to hold, or equivalently, for: 

et<-fe'^^^\ 10M7e'^(^) < e. 

Given e > 0, define 7 so that lOM^e^^^ = e. Then the second inequality 
holds, and so will the first, provided: 

et < 7e-2^-^ = (e/10M)e-^^^, 

which is true for any e > 0, if 5 is defined via 6 := e~^^^/10M (recall 
te[0,S]). 

Note that, since m > at f = 0, it follows that m is negative-definite at 
t = 0, and hence also for small time, and we claim that this persists through- 
out E^ , so that (letting e ^ 0) m < in £'^. Restarting the argument at 
t = 6, we see this is enough to prove the theorem. 
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To prove this claim, suppose (by contradiction) m acquires a null eigen- 
vector / F e R" at a point zi = {xi,ti) G E^, with h e (0,(5] the first 
time this happens. 

Let fiz) := rhijV'V^, z e (that is, we 'extend' V to E^ as a 
constant vector.) It follows from the preceding that / satisfies in E^: 

dtf < trgd^rhijV'V^ + X ■ dm^^Y'V^ + q,^{;m)V'V^ - 

Noting that tVgd'^rhijVW^ = tvgd^f and X ■ rhijV^V^ = df ■ X, and using 
the null eigenvector condition for q, we find that / satisfies in E^ the strict 
inequality: 

dj <trgdy + df-X. 

This shows xi cannot be an interior point of D{ti), for then (as a first-time 
interior maximum point for /) we would have tVgd'^ f{zi) < and df{zi) = 
0, contradicting dtf{zi) > 0. Thus xi G dD{t\). Since / satisfies the 
differential inequality just stated and zi = (xi,ti) is a first-time boundary 
maximum in & , the parabolic Hopf lemma implies dnf{zi) < 0. On the 
other hand, as seen above: 

dnf = dnrhijV'V^ > bij{zi,m{zi))V^V^ > 0, 

from the boundary null-eigenvector condition. This contradiction concludes 
the proof. 

Corollary 13.2. Suppose m G C^'^(£',S) satisfies the same differential 

inequality, with the same hypotheses on the coefficients as in the theorem 
(including the null eigenvector condition for q), and the boundary conditions: 

m{z){n, t) =0, Vz = {x, t) G diE, r G T^dD{t); 

n^nP dnmij{z) = (Vnm)(n,n) > bnniz,m{z)); 

TW^dnniij = (Vnm)(r, r) > 6*»"(z, m(z))(r, r), r G T^dD{t), 

for functions bnn{z,in) from E x S to M. and 6*"" assigning to {z,m), z = 
{x,t), a quadratic form in T^dDit). Suppose 6„„ > in x S and 6*"" 
satisfies: 

rh(r,T) = for some r G T^dD{t) ^ 6*""(z, m)(r, r) > 0. 
Then, as in the theorem, concavity is preserved: 

m<0ati = 0^m<0in^. 
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Proof. This is proved as the theorem, with the following change in the 
last part of the proof: if ^ V £ M" is a null eigenvector of m (defined as 
in the proof of the theorem) at a boundary point zi = (oi, ti) £ diE, write: 

y = y"n + y^, eT^,,dD{ti). 

Assume first V"' ^ 0. Then (noting that m splits at the boundary if 
m does), we see that n is a null eigenvector of m at zi, so we define 
f{z) = mij{z)'n/'{zi)n^ (zi) and repeat the argument. At zi, {Vn'rh){n,n) = 
bnnizi,ifi{zi)) > leads to a contradiction with the parabolic Hopf lemma, 
as before. 

If V"' = 0, then V'^ S TxidD{ti) must be a null eigenvector of m 
at the boundary point Zi, and then we run the argument with f{z) = 
rh{z){V'^ ,V'^), leading to a contradiction, as before. 

Corollary 13.3. For MCM of graphs with constant-angle boundary 
conditions, weak concavity is preserved: 

/i<0att = 0^/i<0in£;. 

Proof. The conditions of the theorem hold, and the expressions obtained 
for Vnh in the preceding section easily imply that the boundary conditions 
in Corollary 13.2 are satisfied; hence the claim follows from Corollary 13.2. 

Remark 13.1. It seems plausible that a slightly different version of the 
result in this section could be used to strengthen the conclusions of [9j . This 
is currently being considered. 

14. An improved continuation criterion. 

In this section we improve the continuation criterion: if sup^ \ h\g = qq 
is finite, the solution can be continued past T. 

The argument given below works in all dimensions, but for simplicity of 
notation we deal here only with the two-dimensional case: St is a surface, 
the moving boundary Tt is a curve in M?. Assuming such a bound on \h\g, 
given the results in section 12 all we have to do is bound {'Vrh){T,T) and 
{Wrh){n, n), where r = is a unit vector field tangent to Tt- At the moving 
boundary: H = (3'^h{n, n)+h{T, r), and we already showed t{H) is bounded, 
so it suffices to bound one of these quantities. 

We adopt the notation: f ^ g ii f — g is bounded in E by constants 
depending only on the initial data and oq. 
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Consider the vector fields in D{t) C M^: 

UJ = —[wi,W2], U = VU'^ = [—W2,Wl]. 

It is easy to verify the following: 

(w, oj)g = 0, \u;\l = \oj\l = \Dw\l := wf + w\. 

Thus we may think of {a;, cD} as a 'conformal pseudo- frame' {uj and uj vanish 
when Dw = 0), defined on all of D($)- Moreover, at the boundary dD{t): 

UJ = Pon, UJ = —n := — r, 

where {t, n} is an euclidean-orthonormal frame along Tf. Thus u and ui 
supply 'canonical' extensions of n, r to the interior of D(t), as uniformly 
bounded vector fields. 

Recall the boundary conditions for h: 

h{LO,ib) = 0, {Vnh){io,io) ^ 0, {Vnh){ib,Lo) ^ on dD{t). 

These give the boundary conditions for the components /in, hi2, /122 of h 
in the standard basis of E?. As shown in appendix 2, these three functions 
are solutions of a linear parabolic system in £^ C x [0, T] (a non-cylindrical 
domain), with bounded coefficients: 

dtihij) - g^\hij)ki + 2h^d^{hkj) + 2h^d^{hik) = Cij. 

The boundary conditions also have bounded coefficients: 

n\''{h22 - hn) + [{n'f - {n''f]hi2 = 0, 

(n^fdnihu) + 2n^n'^dn{hi2) + (n^)^dn(/i22) = h, 

(n^fdnihu) - 2n^n^dn{hi2) + (n^fdn{h22) = 62. 

The only thing left to do is to argue that this set of linear equations and 
boundary conditions define a parabolic system. Then it follows from the 
'global gradient bounds' of linear theory that also the tangential derivatives 
drhij are bounded on diE, which leads quickly to the desired conclusion. 

We need to verify the 'complementarity conditions' hold for this system, 
so we proceed as in Section 6 (up to a point.) Fix a point zq = (yo, ^o) ^ 9iE 
and 'freeze coefficients' there. Consider a manifold-with-boundary chart 
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{y,t) I— (p,(T, s) mapping a neighborhood of zq in E to {p > 0} x R x R+. 
Here p is the coordinate normal to dD(t), a parametrizes dD(t) and slices 
{s = const.} correspond to {t = const.}. 

Let hij{p,a,s) = hij{y,t) be the unknown functions in the new coordi- 
nates. The corresponding system is: 

dghij - {P^(hij)pp + {hij)aa) + 4{hjk)p + Cj{hik)p = tpij. 

Here the cf are constants. The boundary conditions can also be easily 
written down (freeze the n* to their value n^ at zq and replace dn{hij) by 
hp.) It is natural to consider the linear transformation of the unknown 
functions: 

(hll,hi2,h22) ^ (/ll,/l2,/22), fij = fij{p,a,s) 

defined by: 

fii = {nlfhu + 2nlnlhi2 + (^0)^/122 
/22 = {nlfliii - 2nlnlhi2 + (^0)^/122 
/12 = {nlnl){h22 - hi,) + - {nl)']h22. 

Since the principal part of the linear system for the hij is diagonal, the 
principal part of the system for fij is exactly the same (while the lower- 
order terms have different values): 

dsfij - {0^{hj)pp + ifij)aa) + cl{fjk)p + C^jifik)p = Vij, 

This linear transformation is invertible, with inverse given by: 

hn = {nl)^f22 + 2nlnlh2 + (n^)'/ii 

h22 = {nl)^f22 - 2nlnlh2 + {nD^fii 
hi2 = nlnlifu - f22) + [(4)^ - (n^)']/i2 

Thus the original system frozen at zq (for the hij), with its boundary 
conditions, satisfies the complementarity condition if an only if the same 
holds for the fij system, with the transformed boundary conditions. But 
these take a very simple form: 

/l2|p=0 = 0, (^p(/ll)|p=0 = ^1) dp{f22)\p=0 = &2- 

These are standard Dirichlet (resp. Neumann) boundary conditions for a 
standard 3x3 parabolic system (decoupled to highest order). Hence the 
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original system (for the hij) with boundary conditions satisfies 'complemen- 
tarity' at each point of diE. 

In particular, the global gradient estimates hold for the linear system 
(with uniformly bounded coefficients and boundary conditions) in the un- 
knowns hij, and we have: 

\drhij\ < M in diE, 

for any tangential unit vector field r, for some M depending only on gq and 
the initial data. This clearly implies bounds on (\/rh){T, r) and (\/rh){n, n). 
Combining with the results in section 12, we have the following conclusion: 

Proposition 14.1. Assume the maximal existence time Tmax is finite. 
Then: 

limsup sup \h\g = oo. 

t^Tmax dD(t) 

Remark. It should be clear that the argument works in all dimensions; 
this will be included in the final version of the paper. 

15. Final comments. 

1. The main step missing for the global existence result 

lim diamiTit) = 

(in the concave case) is showing that a lower bound on diameter gives an 
upper bound for \h\g. This may follow from properties of the support func- 
tion (based on a point in R" common to all Sj), but remains to be addressed 
(work in progress). If confirmed, this would correspond to Theorem 1 in [5] 
for lens-type curve networks. An issue apparently completely unexplored in 
dimensions above 1 is existence-uniqueness of 'homothetic solutions' for this 
problem. 

2. We state here the local existence theorem for configurations of graphs 
over domains with moving boundaries. In this setting, a triple junction 
configuration consists of three embedded hypersurfaces in M"+-^, 
graphs of functions defined over time-dependent domains D^{t),D'^{t) C 
M" [D^ covered by one graph, by two graphs), satisfying the following 
conditions: (1) The intersect along an (n — l)-dimensional graph A(t) 
(the 'junction'), along which the upward unit normals satisfy the relation: 
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A^i + N2 = N^. (2) If a fixed support hypersurface S C M"^^ is given (also 
a graph, not necessarily connected), the intersect S orthogonally. 

Topologically, in the case of bounded domains one has the following 
examples: (i) ('lens' type) 2 disks (or two annuli) covering -D^(i) and one 
annulus covering D^{t); (ii) ('exterior' type) two annuli covering D'^it) and 
one disk covering D^{t). The boundary component of the annuli disjoint 
from the junction intersects the support hypersurface S orthogonally for 
each t. 

Let (/ = 1,2,3) be graphs of C3+" funct ions over C^^" domains 
Dq,Dq C M"-, defining a triple junction configuration and satisfying the 
compatibility condition for the mean curvatures on the common boundary 
To of and Dg: 

+ = H^. 

Then there exists T > depending only on the initial data, and functions 
e C2+"'i+"/2(Q/)^ c M" X [0,r), so that the graphs of w^{.,t) : 

D\t) R define a triple junction configuration for each t € [0, T), moving 

by mean curvature. 

The proof will be given elsewhere. 

3. An interesting issue we have not addressed here is whether one has 
breakdown of uniqueness for initial data of lower regularity, or if the 'or- 
thogonality condition' at the junction is removed. For curve networks, non- 
uniqueness has been considered in [8]; but neither a drop in regularity (from 
initial data to solution, in Holder spaces) nor the orthogonality condition 
play a role in the case of curves. 

Appendix 1: Proof of lemma 4.1. 

Throughout the proof, n denotes the inner unit normal at dD, extended 
to a tubular neighborhood J\f so that Dnfi = 0. Since D is uniformly C^"*"", if 
follows that n £ C'^"^"(9D), with uniform bounds. Denote by p the distance 
to the boundary (so Dp = n in J\f). Let C ^ C^{D) be a cutoff function, 
with C = 1 in M C AA, C = in Z? \ AA. 

We find (p of the form: 

ip{x) = X + C{x)f{x)n{x) 

with / E C2+"(AA). The 1-jet conditions on Lf at dD translate to the 
conditions on /: 

f\dD = 0, Df\QD = 0, D^f{n, n)\QD = Af^gj^ = h. 
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Now use: 

Lemma A.l. Let be a uniformly (7^+"^ domain with boundary distance 
function p > 0. Let h G C"{dD) be a bounded function. Then there exists 
an extension g G C°°{D) fi C{D) so that g\QD = h, sup^ 1^1 < supg^ \h\ and 

Given this lemma, all we have to do is set / = {l/2)p^g, which clearly 
satisfies all the requirements (in particular, A/ = /i at dD.) 

To verify that is a diffeomorphism, it suffices to check that \C,fn\ci (in 
H d {p < po}) is small if po is small. This is easily seen: 

|C/n|co < (l/2)pg|5lc"; 
\Df\ < (l/2)p^||<7|| 

on J\f , since Df G Ci+"(L») and Df\Qr> = 0. And finally, with A the second 
fundamental form of dD: 

\Dn\ < \A\co => \fDn\ < il/2)pl\g\co\A\co. 



A word about Lemma A.l. (This is probably in the literature, but I 
don't know a reference.) If D is the upper half-space, we solve = in D 
with boundary values h. Then the estimate 

[D\p^P*h)]^''\D)<c\h\c^^on) 

follows by direct computation with the Poisson kernel P; for the rest of the 
norm, use interpolation. Then transfer the estimate to a general domain 
using 'adapted local charts', in which p in D corresponds to the vertical 
coordinate in the upper half-space. (It is easy to see that at each boundary 
point there is a C^"*"" adapted chart, with uniform bounds.) 

Appendix 2: Evolution equations for the second fundamental form. 

We consider mean curvature motion of graphs: 

wt = g^^Wij = vH, V = ^/T+JDw^. 
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In this appendix we include evolution equations for geometric quantities, in 
terms of the operators: 

dt -Ag, L = dt- trgS. 
It is often convenient to use the vector field in D{t): 

uj := -Dw. 

V 

Since —lo is the M" component of the unit normal N and L[N] = \h\'^N, we 
have: 

L[u;'] = \h\y, \h\l:=g"'g^'hijhki. 
Here h = {hij) is the pullback to D{t) of the second fundamental form A: 

h{di, dj) = hij = A{Gi, Gj) = ^Wij. 

First, denoting by V the pullback to D{t) of the induced connection 
(that is, G^{VxY) = V§^^G*y for any vector fields X,Y in D{t)), and 
using the definition: 

Vg^G,- = Gij-{Gij,N)N = [0,Wij]-^Wij[-Dw,l] = ^[Dw,\Dwf] = "^G.Dw, 

we conclude: ^ 

Vo^dj = -hijDw = hijU. 

Prom this one derives easily a useful expression relating the Laplace- 
Beltrami operator and the operator trgdP acting on functions: 

^gf = trgcff Wmfm = tvgd^f - Hdu,f. 



We also have, for the covariant derivatives of h with respect to the eu- 
clidean connection and to V = V^: 

(Here V/i is the symmetric (3, 0)-tensor with components: Vmhij = dm^){di^ ^i)-) 

Iterating this and taking g'-traces yields (using the Codazzi identity and 
the easily verified relation diU^ = /if := g^^hij): 
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+2[hip{h'^)jg + {h%phjg + HhiphjqpPLv'i + 2{h% + 2{h%jh{Lv,Lv). 

Here the powers and of h are the symmetric 2-tensors defined used 
the metric: 

{h%j := g^^hikhpj = h^hpj, {h% := g''^ g^^ hikhpihgj . 
Note also that: 

[h^Vkhjp + h]Vkhip]uP = VUh%j, 
using the Codazzi identity. 

Evolution equations for h. 

Starting from Gt = vHcn+i = H{N + l[Dw, iDwl"^]) = HN + HG^co 
and Nt = -V^H -Hv'^V^v (where V^/ = g'^ fjGi and V/ = g'^ fjdi) we 
have: 

= {{HN)ij,N) - {Gij,V^H) - ^{Gij,V^v) + {{HG,u\j,N). 
Using the easily derived facts: 

{Ni,,N) = ^h\di,dj), 

Hij - {Gij,V''H) = {VdH){di,dj), 
^{Gij,V^v) = h{uj,uj)hij, 

we obtain: 

dtihij) = {VdH){di,dj) - Hh\di,dj) - Hh{u,u)hij + {{HG,io)ij, N), 
where: 

{{HG^u)ij,N) = Hi{{G^u)j,N) + Hj{{G^u)i,N) + H{{G^uj)ij,N). 
To identify the terms, computation shows that: 

{{G,u:)i,N) = h{uj,di), 

and hence, using also: 

V\{G^oj) = G*(Va.a;), Vg.o; = (/if + u}'^hiqUjP)dp = ^ hipdp, 

p 
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we obtain (using aj^dj{hik) = V^hij + 2h{di, u)h{dj , u)): 
{{G,u;)ij,N) = dj{io''hik)-{V^,{G,u;),djN) = h^^hik+^^dj{hik)+h{dj,V a,^) 
= {^ujh)ij + {h^)ij + 2h{u, di)h{u, dj) + ^ hiphj 



VP 

p 

= {y<^h)ij + 2{}?)ij + 3/i(u;, di)h{u, dj), 

since Ylp^ip^^jp = + h{oj,di)h{u!,dj). Combining all the terms yields 

the result: 

dtihij) = {VdH){di,dj) + HV^hij + Hih{cj, dj) + Hjh{u, di) 
+H{h^)ij + 3Hh{u},di)h{u},dj) — Hh{u},u})hij. 

Prom this expression and Simons' identity (in tensorial form): 
VdH = Agh + \h\lh- Hh'^, 
we obtain easily a tensorial 'heat equation' for h: 

[{dt - Ag)h]ij = HV^hij + Hih{io, dj) + Hjhicu, di) 
+\h\ghij + 3Hh{di, ij)h{dj , uj) — Hh{uj,uj)hij. 

Using the earlier computation relating A^/i (the tensorial Laplacian of 
h) and tVgCph, we obtain from this the evolution equation in terms of L: 

L[hij] = -2V^{h%j + C^j, 

Cij := -2[h{di,u)h^{dj,u) + h'^{di,u)h{dj,u)] - 2{h%j - 2{h%jh{u,u) 
+\h^ghij + Hh{di,(jj)h{dj,ijj) — Hh{u,u)hij. 



Time derivatives and evolution equations for lo and g. 

It is sometimes convenient to use the 'Weingarten operator': 

S{X) := S{X'di) = h)X^di. 
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The time derivative of a; is simply minus the time derivative of the M" 
component of N. In addition, one computes easily that ^ = <S'(a;), so we 
have: 

TT 

dtcu = VH + —Vv = VH + HS(io). 

V 

For the metric and 'inverse metric' tensors we have: from dtgij = {wiWj)t 
and Wit = {vH)i: 

dtQij = v'^{HiJ + HjJ) + v'^H{h{uj, di)J + h{uj, dj)J), 

and then, using dtg''^ = -g'^'^dtghig^^: 

dtg^i = -[(yHfJ + {vw)uj'] - H[s{ojyJ + s{ujyJ]. 

Since we know the evolution equation of a;, it is easy to obtain that of 

L[g^^] = -L[a;V] = -L[uj']io^ + 2g^\dkJ){diu^) - JL[u^. 
Using dkUj'^ = Kf., we find: 

L[g^^] = -2\h\liv'J + 2{hy^. 

It is also easy to see that dkg^^ = —{h\u}^ + h-'^oj'^). 

Evolution of mean curvature. 

To compute the evolution equation for H = g^^hij, we just need to 
remember g'^^ is time-dependent: 

{dt-Ag)H = {dtg'^)ihij)+trg[{dt-Ag)h] = -2h{VH,uj)-2Hh\uj,uj)+trg[{dt-Ag)h]. 
The result is: 

{dt - Ag)H = Hd^H + \h\lH + Hh'^iu, u) - H^h{u, u). 

Since L[f] = {dt — Ag)f — Hd^^f (for any /), we see that the equation in 
terms of L has no first-order terms: 

L[H] = \h\]H + Hh^{uj,uj)-H^h{uj,uj) 

Remark. One can also find L[H\ starting from the expression: 

L[g'ihij] = L[g'i]Kj+g'iL[hij] - 2g^\dkg'^){dihij). 

This may be used to check the calculation. 
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Evolution of the Weingarten operator. 

The tensorial Laplacian of S is the (1, 1) tensor AgS with components 
Aghj. We have: 

Agh'; = g'^Aghij, or ((A^^)^,^), = {Agh){X,Y). 

The evolution equation is easily obtained: 

{dt - Ag)hl; = {dtg"')hij + g'^idt - Ag)hij 

= HVu,h'^+Hjh'l'J-Hih^jUj''+\h\lh'^+2HSiio)''h{io,dj)-Hh{uj,uj)h^^^ 
Remark: Since the components of VS" are given by: 

{V^S){dj) = {V^h])dk, V^h] = d^{h]) + h\u, dj)u^ - hiu, dj)S{u)\ 

we see that upon setting j = k and adding over k we recover the evolution 
equation for H. 

The evolution equation for hj in terms of L follows from the calculation: 

L[h'j] = L[g'''\K^+g'^L[h^,] - 25™"(5^<7^')(5„/i.,) 

= -2{V^hl)hJ + {d^h\l)u^ 
+\h\lh] - Hh{u,u)h] + HS{ufh{dj,u) + 2h^{dj,u)u^ - 2{h^)^ujP h{dj , u) . 
Setting j = k and adding over k, we recover the earlier expression for L[H]. 

Evolution of\h\g. 

The fact that g^^ is time-dependent introduces an additional term in the 
usual expression: 

{dt - Ag)\h\l = -2\Vh\l + 2{h, {dt - Ag)h)g + 2 (^^^^ ) (/l^ . 

Using the expressions given earlier, one easily finds: 

{dt - Ag)\h\l = -2\Vh\l + Hd^\h\l + 2|/i|^ - AHh\u^,u) - 2H\h\lh{u,u), 
L[\h?g] = -2|V/i|^ + 2\h\^g - AHh\u,u) - 2H\h\lh{u,u). 



54 



References 

[I] Ecker, K.; Huisken, G. Interior estimates for hypersurfaces moving by 
mean curvature. Invent. Math. 105 (1991), no. 3, 547-569. 

[2] Ecker, K. Regularity theory for mean curvature flow. Progress in Nonlin- 
ear Differential Equations and their Applications, 57. Birkhuser Boston, 
Inc., Boston, MA, 2004. xiv+165 pp. ISBN: 0-8176-3243-3 

[3] Eidelman, S. D.; Zhitarashu, N. V. Parabolic boundary value problems. 
Operator Theory: Advances and Applications, 101. Birkhuser Verlag, 
Basel, 1998. xii+298 pp. ISBN: 3-7643-2972-6 

[4] Guan, B. Mean curvature motion of nonparametric hypersurfaces with 
contact angle condition. Elliptic and parabolic methods in geometry 
(Minneapolis, MN, 1994), 47-56, A K Peters, Wellesley, MA, 1996. 

[5] Evolution of convex lens-shaped networks under curve shortening flow 
O. Schniirer, A. Azouani, M. Georgi, J. Hell, N. Jangle, A. Koeller, T. 
Marxen, S. Ritthaler, M. Saez, F. Schulze, B. Smith, (Lens Seminar, 
FU Berlin 2007), arXiv:0711. 11081 

[6] Baconneau, O.; Lunardi, A. Smooth solutions to a class of free bound- 
ary parabolic problems. Trans. Amer. Math. Soc. 356 (2004), no. 3, 
987-1005. 

[7] Mantegazza, C., Novaga, M., Tortorelli, V., Motion by curvature of 
planar networks. Ann. Sc. Norm. Super. Pisa CI. Sci. 5(3) (2004) no. 
2, 235-324. 

[8] Mazzeo, R., Saez, M. Self-similar expanding solutions of the planar 
network flow, ar Xi v : 704 ^ 3 11 3 

[9] Stahl, A. Convergence of solutions to the mean curvature flow with a 
Neumann boundary condition. Calc. Var. Partial Differential Equations 
4 (1996), no. 5, 421-441. 

[10] Solonnikov, V. A. Lectures on evolution free boundary problems: clas- 
sical solutions. Mathematical aspects of evolving interfaces (Funchal, 
2000), 123-175, Lecture Notes in Math., 1812, Springer, Berlin, 2003. 

[II] Struwe, M. The existence of surfaces of constant mean curvature with 
free boundaries. Acta Math. 160 (1988), no. 1-2, 19-64. 



55 



